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Azimuthal angle and rapidity dependence of di-hadron correlations in QCD 



I discuss novel QCD phenomena recently observed in p+p, p+A and A+A collisions, that result from the non- 
linear dynamics of small-x gluons. I focus on di-hadron correlation measurements, as opposed to single-hadron 
observables often too inclusive to distinguish possible new effects from established mechanisms. Specifically, I discuss 
angular correlations of forward di-hadrons in d-\-Au collisions and long-range rapidity correlations in high-multiplicity 
p-hp and Au-hAu collisions. 
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1. Introduction 

The large amount of experimental data collected at RHIC during the last decade, and at the LHC during the last 
year, allows to explore novel QCD eff'ects. In particular, RHIC measurements present a number of features which can 
be consistently interpreted as consequences of a large gluon density in the gold nucleus, and of the corresponding non- 
linear QCD dynamics. Manifestations of this so-called saturation regime of QCD include the suppression of forward 
particle production 1 1 , 2| and the azimuthal de-correlation of forward hadron pairs | 3, 4 | in d-hAu compared to p-hp 
collisions. The recent start of the LHC significantly enhances the possibilities to explore the saturation regime of QCD 
already in p-hp collisions, as shown by the recent observation of long-range rapidity correlations in high-multiplicity 
events |5|. 

We argue below that these data are well described by the eff'ective theory called the Color Glass Condensate (CGC), 
the best approximation of QCD devised so far to quantify parton saturation (see e.g. |6|). It properly takes into ac- 
count both the non-linear evolution of hadronic/nuclear wave functions, and multiple parton interactions in scattering 
processes. Such eff'ects are relevant when the gluon densities involved are large enough, corresponding to partons 
with a small-enough energy fraction x, and the CGC describes reliably only that small-x part of the hadronic/nuclear 
wave function. And although the CGC describes non-linear QCD dynamics, it is still a weakly-coupled theory. 

In p-hA collisions, the saturation regime is better studied with forward particle production, sensitive only to large- 
X partons on the proton side, but mainly to small-x gluons on the nucleus side. Since the large-x part of a proton 
wave function is well understood in perturbative QCD, this allows to investigate the small-x part of the nucleus 
wave function. In A-hA and high-multiplicity p-i-p collisions, the dynamics of small-x gluons can be studied using 
long-range rapidity correlations, sensitive to the early times after the collision rather than the following complicated 
space-time evolution of the system. While the existence of the QCD saturation regime is, at a theoretical level, clear, 
the real challenge from a phenomenological point of view has been to assess to what extent it is relevant at present 
energies. Such was a difficult task, since diff'erent aspect of the parton dynamics may concur in data, and also because 
the limit of asymptotically high energy in which the CGC formalism is developed is not fully realized in current 
experiments. 
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Nevertheless, this has been achieved with the successful simultaneous description of the suppression of particle 
production |7 1 and azimuthal correlations 1 8 1 at forward rapidities in d+Au compared to p+p collisions, using the 
most up-to-date theoretical tools available in the CGC approach. This is explained in Section [2j in parallel with more 
standard mechanisms that successfully describe particle production at mid-rapidity. In Section |3j we discuss long- 
range rapidity correlations in high-multiplicity p-hp and Au-hAu collisions, and argue that the angular structure seen 
in the data and called the ridge may be interpreted as a consequence of large parton densities ||9l[T0l, albeit in these 
cases several other possible interpretations exist. 



2. Di-hadron correlations, p+p vs d+Au collisions 



In the case of double-inclusive hadron production pA hih2X, the partons in the proton (nucleus) wave function 
that can contribute to the cross section carry a fraction of longitudinal momentum bounded from below by Xp (xa), 
with 



\Pi^\ey^^\P2^\ey 



\pii.\e-y' +\p2^\e- 



Xp = ^= , Xa = ^= , (1) 



/snn ySNN 
where /^i^, P2± and ji, y2 denote the transverse momenta and rapidities of the final-state particles, and ^/s^ the 
collision energy per nucleon. One can distinguish three interesting situations, where diff'erent QCD dynamics is at 
play. First, the production of high-pr di-jets, initiated by large-x partons and kinematically bound to be produced at 
mid-rapidity. Second, the production of low-pr di-hadrons at mid-rapidity, sensitive to power corrections and values 
of X still moderately large. Note that moving one of the hadrons to forward rapidities (central-forward case) increases 
significantly the value of Xp compared to the central-central case (for which Xp ^ xa - PtI ^sI^nn)^ but decreases xa 
only marginally; for this reason we shall refer to these two situations as mid-rapidity ones. Finally, the production 
of di-hadrons at forward rapidities (forward-forward case) is sensitive to power corrections in the small-x regime, i.e. 
to saturation. Only the forward-forward case is sensitive to values of xa as small as in the single-inclusive forward 
hadron production case: Xp<\ and ^ 1. 



2.7. High-pr di-jets and the leading-twist approximation 

Fig. [T] shows the cross section for high-;?^ di-jets in p-Fp colli- 
sions at the Tevatron, as the function of the azimuthal angle between 
the jets A0 |11|. The leading-twist approximation of QCD, in which 
the cross section is computed to leading-power of the hard scale, de- 
scribes the data very well, and in this regime the di-jet azimuthal 
correlation is extremely peaked around A0 = tt. It is often assumed 
that this coUinear factorization approach also holds in p-^A collisions, 
meaning that highly-virtual partons in nuclei behave independently as 
they do in protons. Whether or not this is the case could be tested at 
the LHC, starting with inclusive-jet/hadron measurements |[T2l . 

2.2. Low-pr di-hadrons and power corrections at large x 

We shall first focus on di-hadron production at mid-rapidity, in- 
cluding both the central-central and central-forward situations, and 
study the A0 dependence of the double-inclusive hadron production 
cross section, where Acp is the diff'erence between the azimuthal an- 
gles of the measured particles hi and h2. Nuclear eff'ects on di-hadron 
correlations are typically evaluated in terms of the coincidence prob- 
ability to, given a trigger particle in a certain momentum range, pro- 
duce an associated particle in another momentum range. 

In a p-hp or p-hA collision, the coincidence probability is given by 

CP(A(f>) = Npair(A(p)/Ntrig, with 
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Figure 1: High-pj mid-rapidity di-jet cross-section 
measured by the DO collaboration fTTj as the func- 
tion of A0. Although decreasing the pj range 
leads to azimuthal decorrelation, the signal remains 
sharply peaked around A0 = n. This is in agreement 
with leading-twist NLOQCD predictions, which con- 
tain up to four partons in the final state. 



C. Marquet/ Nuclear Physics A 00 (2011) 7-[S| 



3 




Figure 2: The coincidence probability at mid-rapidity as a function of A0. RHIC data show that in the central-central case the away-side peak is 
similar in d-i-Au and p-i-p collisions. In the central-forward case, the Glauber-eikonal (left plot from 1 14 |) and CGC (center plot from 1 15 1) calcu- 
lations predict that the away-side peak is suppressed in d-i-Au compared to p-i-p collisions. This is in agreement with central-forward preliminary 
IdAu data (right plot from 1 16 1), which show that in central d-i-Au collisions, the integral of the coincidence probability is about half that in p-i-p 
collisions, reflecting the depletion of the away-side peak. CoUinear factorization cannot reproduce this behavior, while both Glauber-eikonal and 
CGC calculations predicted it. 



First measurements were performed at RHIC by the PHENIX and STAR collaborations f?,T3l. In the central-central 
case, the coincidence probability features a near-side peak around A0 = 0, when both measured particles belong to 
the same mini-jet, and an away-side peak around A0 = tt, corresponding to hadrons produced back- to-back. In the 
central-forward case, there is naturally no near-side peak. Either in p-hp or d-i-Au collisions, the sizable width of the 
away- side peak cannot be described within the leading-twist collinear factorization framework. This indicates that 
power corrections are important when \pj_\ ~ 2 GeV. At such low transverse momenta, collinear factorization does not 
provide a global picture of particle production at RHIC, even at mid-rapidity. 

There are two main formalisms to take into account power corrections to the leading-twist approximation, the 
Glauber-eikonal approach, and the CGC framework. The former relies on a resummation of incoherent multiple scat- 
terings, and therefore is not applicable at small x, where multiple scatterings are coherent. Performing the complete 
resummation including energy-momentum conservation is a challenging task, and more often than not a detour to the 
strict calculation is taken by introducing non-perturbative parameters such as the so-called intrinsic transverse mo- 
mentum. By contrast, the CGC is valid only in the small-x limit. Coherent multiple scatterings are taken into account 
along with the equally important non-linear parton evolution, since both parametrically contribute to the same level 
when a large gluon density is reached, and including one without the other is not consistent. In both frameworks, 
the multiple scatterings induce a /7r -broadening, characterized either by a non-perturbative scale at large x, or by the 
saturation scale at small x, as explained in the next section. 

Both the Glauber-eikonal 1 14] and CGC 1 15 1 approaches qualitatively describe the di-hadron azimuthal correlation 
data at mid-rapidity, including the broad width of the away- side peak and its depletion in d-hAu collisions when going 
from central-central to central-forward production. This is illustrated in Fig.|2j left and center plots. Such a depletion 
does not occur in p-hp collisions, it is due to nuclear-enhanced power corrections, and the p-hA to p-Fp ratio of the 
integrated coincidence probabilities 

/ JA0 CPpA(A(p) 

tpA = -T (3) 

/ dAcp CPpp(A(P) 

is therefore below unity. In Fig.[2j right plot, recent PHENIX data on I^au l[T6ll are displayed as a function of centrality. 

At the moment, it is not clear whether the physical origin of the /7r -broadening and of the azimuthal de-correlation 
is parton saturation rather than incoherent multiple scatterings. This is probably due to the kinematic window probed 
in these measurements: for a hadron momentum of \p_i_\ ~ 2 GeV, one is sensitive to xa ~ 0.01 to 0.1 at RHIC 
energies. In this regime different physical mechanisms may concur, and neither underlying assumption of the two 
different descriptions is completely fulfilled. The approximations made in CGC calculations will apply best in the 
forward- forward case discussed next. 



C. Marquet/ Nuclear Physics A 00 (2011) 7-[S| 



4 



p+p 7fV+X, Vs = 200 GeV 
r Pt,ll>2 GaV/c, 1 GeV/c<p„<pT,L 
L <7]^>=3.2, <%>=3,1 




d+Au 
0.03 




> ttV+X, Vs = 200 GeV, 2000<IQbbc<40 
Pt,l>2 GeV/c, 1 GeV/c<pT,s<PT,L 
I, <7]l> = 3.1, <?is>=3.2 

m 



Peoks 

0.44±0.02 
n 1,63±0.29 



Preliminary 







p^^>2GeV/c 
: 1 GeV/c < Pt 3<Pt L 


STAR PRELIMINARY 
. p+p (-0.0045) 
. d+Au central (-0.0145) 


« 

• 


\*'. * 


- M • ■ M 









A(|) 



Figure 3: The coincidence probability as a function of A0 for p+p (left) and central d+Au (center) collisions. These preliminary data from 1 3 1 
show a striking nuclear modification of di-hadron azimuthal correlations. The away-side peak, corresponding to hadrons emitted back-to-back, 
is prominent in p+p collisions but is absent in the central d+Au case. Such production of mono-jets was anticipated in the CGC as a signal of 
parton saturation, the data are compared with CGC predictions on the right plot; the away-side peak in p+p is qualitatively described by the CGC 
calculation while its disappearance in central d+Au is quantitatively consistent with the prediction. 



23. Forward di-hadrons and parton saturation at small x 

At RHIC mid-rapidity observables are contaminated too much by large-x^ physics to quantitatively explore sat- 
uration effects. As explained in the introduction, particle production in the forward-rapidity region offers a cleaner 
opportunity. This was realized in ifTVl . where the importance of the forward di-hadron measurement was emphasized, 
and predictions of the azimuthal angle distribution were made. Although these results were only qualitative since for 
instance parton fragmentation was not included, they have been confirmed by subsequent RHIC measurements. The 
recent data collected by the STAR collaboration for the coincidence probability obtained with two neutral pions are 
displayed in Fig.[3j for both p-hp (left plot) and central d-hAu (center plot) collisions. The nuclear modification of the 
di-pion azimuthal correlation is quite impressive, the prominent away-side peak in p-hp collisions is absent in central 
d-hAu collisions, in agreement with the behavior predicted in 1 17 |. 

Let us outline this calculation. As external probes propagate through a target hadron or nucleus, their multiple 
scatterings off the small-x gluons is coherent because the wavelength of these gluons increases with decreasing x, 
eventually becoming bigger than the longitudinal extent of the target. The momentum scale below which such co- 
herent multiple scattering are important is the saturation scale 2^(x), the same scale which characterizes the onset of 
non-linear effects in the hadronic/nuclear wave function. Both are taken into account in the CGC approach, which 
resums all power corrections that are dominant in the small-x limit. Non-linear parton evolution generates an intrinsic 
transverse momentum while some is also gained through multiple scatterings. The magnitude of both is controlled by 
Qs, which is an increasing function with decreasing x. 

In the leading-log(l/x) approximation, the evolution of the saturation scale with x, and more generally of the CGC 
wave function, can obtained from a functional renormalization group equation. In the large-A^c limit, it reduces to the 
non-linear Balitsky-Kovchegov (BK) equation ifTSl [T9l . The recent determination of running coupling corrections 
to the original leading-logarithmic equation ||20l[2Tl has proven an essential step in promoting the BK equation to a 
phenomenological tool. Indeed, the derivation of the rcBK equation, along with numerical simulations I'S?], made 
robust quantitative calculations possible. For instance, after the two parameters of the theory were constrained by 
single-inclusive forward hadron production Q (they are xq, the value of xa below which one starts to trust, and 
therefore use, the CGC framework, and the value of the saturation scale at the starting point of the evolution 2^(xo)), 
parameter- free predictions for the coincidence probability CP(A(p) could be made |8 1. 

In Fig. [3] right plot, these predictions are compared with the STAR data. We see that the disappearance of the 
away-side peak in central d-hAu collisions, compared to p-hp collisions, is quantitatively consistent with the CGC 
calculations. These are only robust in the d-hAu case (for which 2^(xo = 0.02) = 0.6 GeV^), but the extrapolation to 
the p-Fp case (for which Qj(xo = O.Ol) = 0.2 GeV^ is uncomfortably close to A^^^) is displayed in order to show that 
it is qualitatively consistent with the presence an the away- side peak. Since uncorrelated background has not been 
extracted from the data, the overall normalization of the data points has been adjusted by subtracting a constant shift, 
as indicated on the figure. 
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We recall that di-hadron correlations at mid-rapidity, which are sensitive to larger values of xa, feature an away- 
side peak in both p-hp and central d-hAu collisions. The fact that the away- side peak disappears in the central d-hAu 
case, when going from central to forward rapidities, indicates that the effect is correlated with the growth of the nuclear 
gluon density with decreasing xa- Similarly, we predict that in the forward-forward case, the away-side peak reappears 
when going from central to peripheral collisions, or for higher transverse momenta of the measured particles. We are 
not aware of any descriptions of this phenomena that do not invoke saturation effects. We note that another successful 
description based on the KLN saturation model was also proposed 1231 . There, although different assumptions are 
used, the existence of the saturation scale is the crucial ingredient to successfully reproduce the data. 

While more differential measurements of the coincidence probability, as a function of transverse momentum, 
rapidity or centrality, will provide further quantitative tests of the CGC, this piece of evidence strongly indicates 
that we have observed manifestations of the saturation regime of QCD at RHIC. However, there are still progresses 
to be made on the theoretical side. For instance, an approximation was used in [17^ in order to express the 
double-inclusive hadron production cross section (also derived in ||24l [25] |26l and in 1271 with in addition gluon- 
initiated di-jets) in terms of only the single-gluon distribution function, already well constrained. The validity of this 
approximation has been critically examined [28], as it does not allow to correctly implement the non-linear QCD 
evolution, even in the large Nc limit. Still, in practice, no better solution has been proposed so far. This will require 
to numerically solve the x evolution of multi-gluon distributions, which implies to introduce additional parameters 
to characterize the corresponding initial conditions. It was pointed out recently that in the limit of small di-hadron 
momentum imbalance, only single-gluon distributions are needed |29|, but in general this is not the case. 

3. Long-range rapidity correlations, p+p vs Au+Au collisions 

Let us now discuss high-multiplicity p-hp and heavy-ion collisions. When considering low-pr particle production, 
only the small-x part of the wave functions of the colliding objects is relevant, and therefore such collisions can be 
described by the scattering of two CGCs. The system created immediately after such a collision is called the Glasma, 
it provides a weak-coupling description of the early stages after a high-multiplicity p-hp or A-hA collision, before 
evolving into a thermalized quark-gluon plasma it the system size is big enough, as is the case for not-too-peripheral 
heavy-ion collisions. 

If we denote pi and p2 the strong color charges of the incoming hadrons/nuclei, the color field describing the 
dynamics of the small-x gluons is solution of the Yang-Mills equations 

[£)^,Fn=<5>i+<5->2. (4) 

The Glasma field, after the collision, is non-trivial (30): it has a strong component (A^ ~ 1/^^), a component which 
is particle like (A^ ~ 1), and components of any strength in between. To understand how this pre-equilibrium system 
thermalizes, one needs to understand how the Glasma field decays into particles. Right after the collision, the strong 
field component contains all modes. Then, as the field decays, modes with pr > I It are not part of the strong 
component anymore, and for those a particle description becomes more appropriate. After a time of order 1/2^, 
this picture breaks down, and it has been a formidable challenge to determine weather a fast thermalization can be 
achieved within this framework. 

A problem which can be addressed more easily is particle production. As discussed below, a factorization frame- 
work has been established to perform such calculations 1311 . Before we outline it, let us point out that two-particle 
correlations are especially interesting, because they can probe the dynamics of the early times after the collision, if the 
rapidity difference between the two final-state particles is large enough |32|. Indeed, causality imposes that particles 
separated by a rapidity At; can be correlated only by events which happened at 

T < Tf.o. e-^^l^ , (5) 

where the freeze-out time t/.o. denotes the time of last interaction. This is explained in Fig. |4] left diagram, which 
shows the space-time picture of the collision: in the forward light-cone, lines of constant proper time r = ^|x^x~ 
are hyperbolae and lines of constant rapidity ^ = | log(x^/x") are straight lines from the origin. Therefore, while 
heavy-ion collisions are in general not great probes of parton saturation due to the complicated space-time evolution 
of the system, long-range di-hadron correlations are particularly well-suited observables for that purpose. 
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Figure 4: Left: the space-time location of events that may correlate two particles is the intersection of their past Hght-cones. Correlations between 
particles widely separated in rapidity are due to early-time dynamics. Right: typical leading-order diagram for particle production in the Glasma, 
multiple partonic interactions are crucial when low values of x are being probed in the nuclear wave functions. 



3.1. Particle production in the Glasma 

To compute particle production in the collision of two CGCs, the first task is to express the cross sections in terms 
of the Glasma field A^[pi,p2], taking into account multiple partonic interactions, as pictured in Fig.[4j right diagram. 
Then, one needs to average the cross sections over all possible field configurations, distributed according to functional 
probabilities that we denote Wjc^[pi] and Wjc£[p2]- 

{0} = J Dp,Dp2W,MWxM]0[A^]. (6) 

This factorization has been recently proven [ 31 1 to leading order in as, and to all orders in gsPi and gsPi- The small-x 
evolution is contained in the CGC "wave functions" Wx, while the momenta of the final-state particles enter in the 
observables O. For instance, in the case of single-inclusive gluon production, the cross section reads 



d^p^dy 



An aspect of hadronic collisions which happens to be much simpler when large parton densities are involved, com- 
pared to small ones, is multi-gluon production. Indeed, in a collision of dilute objects of small partonic densities, the 
produced gluons predominantly all come from the same ladder, because this is what involves the parton distributions a 
minimum number of times (once for each colliding hadron). By contrast, in a collision of dense objects, the dominant 
contribution comes about when the produced gluons all originate from diff'erent ladders, so as to involve both gluon 
distributions a maximum number of times (a pair for each produced gluon). The multi-gluon production cross section 
then reads: 

dN dN dN 

d^pi . . . d^pn d^pi d^pn 

This power counting argument that qualitatively distinguishes multi-particle production in high-multiplicity p-Fp and 
A-hA collisions compared to minimum-bias p-i-p collisions, also applies when ISj] is large. 

In the leading-log(l/x) approximation, the evolution of W^lp] with decreasing x is given by a functional renor- 
malization group equation. As both the CGC wave function and the Glasma field A^ first have to be evaluated 
numerically, performing the p integrations while fully implementing the non-linear QCD evolution has proven to be an 
extremely challenging task, and has not been achieved yet. In phenomenological studies of two-particle correlations, 
it is instead assumed that W^lp] is a Gaussian distribution, whose variance is evolving with x such that the evolution 
of the single-gluon distribution is correctly reproduced (but that of multi-gluon distributions is not). In this case, the 
double-inclusive hadron production cross section can be expressed in terms of only single-gluon distributions lO, al- 
lowing a practical phenomenological analysis. However, similarly to what was mentioned at the end of Section[2j this 
approximation does not correctly implement the non-linear QCD evolution, as it misses a large number of leading-A^^ 
contributions |[33l . It was also argued recently that, as a consequence of using Gaussian distributions, two-particle 
correlations are underestimated 041 . 
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Figure 5: Left: figure showing that particles emitted from the same flux tube are correlated in rapidity in a nearly boost-invariant way, and explaining 
how radial flow turn this a priori A0-independent correlation into a ridge structure. Center: STAR measurement of di-hadron correlations in central 
Au+Au collisions 1 37 1, featuring long-range rapidity correlations on both sides of the jet peak, in qualitative agreement with the CGC expectations. 
Right: comparison between PHOBOS data 1 36 1 and GCG calculations | 9 1, of the magnitude of the ridge at A0 = as a function of A?/. 

3.2. The ridge in p-\-p vs Au+Au collisions 

Two of the main features of the Glasma field are that, at leading-order it is //-independent, and that its correlation 
length in the transverse plane IjQs is much smaller than hadronic/nuclear radii. Consequently, the color field A^ 
is made of boost-invariant longitudinal flux tubes of transverse size 1/2^. Because of this flux-tube structure, the 
Glasma is a natural candidate to explain the ridge-shaped two-particle correlations observed in Au-hAu collisions at 
RHIC EH EH [521. This is illustrated in Fig. [5] left drawing: the particles can come from difl'erent tubes in which case 
they are uncorrelated, or come from the same tube which creates an //-independent correlation. 

Due to the strong radial flow induced at a later stage by the collective motion of the system, this correlation is 
turned into a ridge collimated around A0 = 0, such as the one displayed in Fig. [5] center plot, obtained from a di- 
hadron correlation measurement in central Au-hAu collisions. In agreement with our interpretation, this ridge-shaped 
structure is only observed with particles of low-enough transverse momenta, that participate in the collection motion. 
It disappears in peripheral collisions as well, since these do not create a quark-gluon plasma, and therefore any radial 
flow. In this case, only a strong-enough intrinsic A0 dependence of the two-particle correlations could reveal the 
long-range rapidity correlations. 

The right plot of Fig. [5] shows the magnitude of the ridge as a function of A//, for A0 = 0. The data are successfully 
compared with a calculation that combines an estimate of the jet peak around A// = 0, with the nearly //-independent 
correlation computed from the double-inclusive particle production cross section in the CGC |9 |. The figure also 
shows that the ridge is several units long in rapidity, confirming that these correlations originate at early times. 

Using high-multiplicity events in p-hp collisions at the LHC, the CMS collaboration |5 1 recently answered the 
question whether, in a system without much collective flow, one could detect the rather- small intrinsic A0 dependence 
of di-hadron correlations in QCD at high energies. It happens to be in qualitative agreement with CGC expectations 
ifTOI. For instance, the left diagram of Fig. [6j which involves two (A,B) pairs of gluon distributions as previously 
explained, is responsible for a A0 dependence of two-particle correlations, when the transverse momentum of both 
final-state hadrons is of order Qs. This QCD correlation is irrelevant in heavy-ion collisions, where it is overwhelmed 
by radial flow, however it dominates when the system is too small to develop any collective motion. This seems to 
be the case in high-multiplicity p-hp collisions at the LHC, where no ridge is seen for XosN-pr particles (Fig.[6j center 
plot). On the contrary, the ridge is indeed present when pr ~ Qs (Fig.[6j right plot). 

4. Conclusions 

To conclude, let us emphasize that saturation-based approaches were the only ones to correctly predict the sup- 
pression of forward particle production [ 38 , 39] and the azimuthal de-correlation of forward hadron pairs 1 17 1 in p-hA 
vs p-hp collisions at RHIC, as well as the ridge phenomenon in high-multiplicity p-hp collisions at the LHC BOl . 
While alternative, subsequent-to-data and distinctive explanations of these difl'erent phenomena have been and may 
still be proposed, we are not aware of any formalism that can describe all of them in a global way, other than the CGC, 
which, let us not forget, is also known to describe well all small-x sensitive observables in e-hp collisions at HERA. 
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Figure 6: Left: a typical diagram that generates an angular collimation around A0 = nearly independent of Atj, the rapidity separation between 
the two hadrons. However, it is non-zero only for hadron transverse momenta comparable to Qs, and it is dominant only at large A77 away from the 
jet peak. Center and right: CMS measurement of di-hadron correlations in high-multiplicity p-i-p collisions | 5 1, as a function of A0 and A77. The 
data is consistent with the CGC picture, as the so-called near-side ridge is only observed when 1 < pr < 3 GeV. 
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